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ABSTRACT 

The relative importance of different initiation mechanisms for coronal mass ejections 
(CMEs) on the Sun is uncertain. One possible mechanism is the loss of equilibrium of 
coronal magnetic flux ropes formed gradually by large-scale surface motions. In this 
paper, the locations of flux rope ejections in a recently-developed quasi-static global 
evolution model are compared with observed CME source locations over a 4.5-month 
period in 1999. Using EUV data, the low-coronal source locations are determined 
unambiguously for 98 out of 330 CMEs. An alternative method of determining the 
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source locations using recorded Ha events was found to be too inaccurate. Despite 
the incomplete observations, positive correlation (with coefficient up to 0.49) is found 
between the distributions of observed and simulated ejections, but only when binned 
into periods of one month or longer. This binning timescale corresponds to the time 
interval at which magnetogram data are assimilated into the coronal simulations, and 
the correlation arises primarily from the large-scale surface magnetic field distribution; 
only a weak dependence is found on the magnetic helicity imparted to the emerging 
active regions. The simulations are limited in two main ways: they produce fewer 
ejections, and they do not reproduce the strong clustering of observed CME sources 
into active regions. Due to this clustering, the horizontal gradient of radial photospheric 
magnetic field is better correlated with the observed CME source distribution (coefficient 
0.67). Our results suggest that, while the gradual formation of magnetic flux ropes over 
weeks can account for many observed CMEs, especially at higher latitudes, there exists 
a second class of CMEs (at least half) for which dynamic active region flux emergence 
on shorter timescales must be the dominant factor. Improving our understanding of 
CME initiation in future will require both more comprehensive observations of CME 
source regions and more detailed magnetic field simulations. 

Subject headings: Sun: corona — Sun: coronal mass ejections — Sun: evolution — 
Sun: magnetic fields 



1. Introduction 



Coronal mass ejections (CMEs) are a major a rea of interest in so l ar physics, bot h because 



of their influence on the near-Earth environment (iGosling et al 



cause of their role in the global magnetic field evolution ( Bieber &: Rust 19951 : Zhang &: Low 2001 



1974 



Schwenn 



200i 



and b 



Owens et al.ll2007l ). However, the physical mechanism or mechanisms responsible for their initiation 
remains an open question; while there is agreement that the driving energy must originate in the 
magneti c field, the manner in which this free ma gnetic energy is built up and released is still under 



debate (|Forbesll200d : lKhmchukll200ll : ILowI 12001 



The present study considers a particular model for CME initiation: the quasi-static build- 
up and sudden loss of equilibrium of coronal magnetic flux ropes, in response to the large-scale 
surface motions of differential rotation and meridional flow, and to flux cancellation. These 
motions generate electric currents in the corona and the resulting magnetic helicity becomes 
conce ntrated iri twis t ed flux rope structures above po larity inversion lines in the photospheric 
field (|Pneumanlll983l : Ivan Ballegooijen &: Martenall989l ). This formation process takes place over 
timescales on the order of a month and is therefore an appealing mechanism for the formation of qui- 



escent filaments (Ivan Balleeooiien &: Martens 



1990 



Mackay &: van Ballegooijen l200,4 ICibson fc Fanlbood ). Recent simulations of the global corona 



Zirker et al. 



1997 



van Ballegooijen et al. 



199S 



- 3 - 



demon strate that this mechanism can account for th e observed hemispheric pattern of filament chi- 
rahtv (lYeates et al.ll2008al : lYeates &: Mackavll2009al ). as weU as leading to the ejection of flux ropes 
([Yeates &: Mackavll2009bl ). On shorter timescales (hours to days), observations show that magnetic 
flux in active regi ons can emerge highly twiste d, carrying substantial currents from beneath the 
photosphere (e.g., iLeka et al.l Il996l : lLiteal2009l ). Although we include a net magnetic helicity in 
emerging active regions, our large-scale, quasi-static model cannot follow the dynamic evolution 
of the flux emergence and its rapid reconfiguration on entering the corona. Indeed, the model 
was originally designed to follow the large-scale dispersal of magnetic flux and helicity from active 
regions, so does not, at present, follow the detailed evolution on short timescales that is observed 
inside active regions. This paper aims to determine how the distribution of flux rope ejections in 
the model compares wit h that of observed CMEs. We do this through a direct comparison between 
the simulation results of lYeates &: Mackavl (j2009bl ) and CME observations over a 4.5-month period 
in 1999 using the Solar & Heliospheric Observatory (SOHO). This is the first such comparison of 
CMEs with long-term simulations representing observed magnetic configurations on the Sun. 

A difficulty with such a direct observational comparison lies in the identification of CME source 
regions in the low corona. There are multiple low-coronal signatures that may indicate the CME 
source region, including: coronal dimming regions, a post-eruptive arcade, a filament eruption 
and/or a coronal wave. However , on occasion, CMEs ar e detected which appear to have no low- 
coronal source region at all (e.g. iRobbrecht et al.ll2009bl ). Such CMEs may originate from higher 
in the corona, and lack a low-coronal signature. Even when a low-coronal signature is observed, 
confidently linking it with a particular CME requires a case-by-case study, and consideration of both 
spatial location and a plausible temporal association between the CME and its low coronal signature 
is required. Part of our aim in this paper is therefore to consider how meaningful a comparison of 
simulations with observations can be, and what would be needed for a better comparison. 

The key feature of this work is that we simulate the global magnetic field in the solar corona, 
so that flux ropes form in a time-dependent manner at different locations on the Sun, in response 
to emergence of active regions and large-scale surface motions. This is in contrast to previous 
studies which have typically modelled a single CME event in a simplified magnetic configuration, 
in order to consider the basic physical processes leading to loss of equilibrium. These studies 
indicate that the evolution of flux rop es, and whether sudden eruption will occur, depends both 
on the photqspheric footpoint motions (IForbes &: Isenbergjll99ll : iMikic &: Linkeiill994l : lAmari et al 



1996 



2003 



(jlsenberg et al. 



1993 



Mackav &: van BallegooiienI bood) and on the overlying, bac kground magnetic field 



Antiochos et al.lll999l : iLin &: van BallegooiienI |2005| ). In the model used in 



this paper, the footpoint motions are determined by the large-scale surface motions of differential 
rotation, meridional fiow, and supergranular diffusion, while the overlying magnetic field at the 
location of each flux rope is determined self-consistently in the global magnetic conflguration. This 
allows us to place constraints on the applicability of this CME initiation mechanism on the real 
Sun. In this context, several recent studies have looked at the large-scale magnetic topology of 
observed CME source regions, focusing on whether the background field is bipolar or quadrupolar, 
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because CME initiation model s differ fundament ally in this respect (jLi &: Luhmannll2006l : iBarnes 
20071 : lUgarte-Urra et alj 120071 : ICook et al.l l2009l ). Our model takes into account this topology 
automatically, because both the flux rope and its overlying field are part of the global magnetic 
configuration. 

The lack of any previous such global models may be explai ned by the need for c oronal electric 
currents, required in order to store free magnetic energy (e.g., ISchrijver et al.l 120081 ). These must 
be built up in a time-dependent manner, either rapidly through the emergence of pre-twisted struc- 
tures, or over longer timescales by surface shearing, as in the simulations in this paper. Global 
full-MHD models are too computationally demanding to model this time evo lution over many 



weeks, although they have been succes sful in modelling both global equilibria (IRilev et al 



Cohen et al 



2008 



2007: 



Lionello et al. 



20091') and the evolution of individual CME events (jRiley et al 



2006 



sT Manchester et al.., 20081 : Cohen et al. 2009 ). Instead, our model uses a quasi-static approxi 



mation to the magnetic field evolution in order to follow the formation of currents and the transport 
of magnetic helicity, albeit in a simplified manner. 

In Section [2] we describe the instrumentation and data reduction used to identify the source 
regions of observed CMEs. The main features of the simulations are outlined in Section [3l before 
comparing the observed and simulated distributions of CME sources in Section [H The relation 
between the two is discussed in Section [5l and conclusions in Section [H where we also make 
recommendations for future study. 



2. Observations of CME Source Regions 



We began by compiling a list of CMEs between 199 9 May 13 and 1999 September 26 from the 
CDAW (coordinated data analysis workshops) catalog ( Yashiro et al. 2004). This is the standard 
manually-compiled list of CMEs observed by the Larg e Angle and Spectrometric Coronagraph 
(LASCO) experiment since 1996 (jBrueckner et al.lll995l ). The C2 and C3 coronagraphs observe 
the white-light corona from 2Rq to 7Rq and 3.7Rq to 32i?0 respectively. Our comparison period 
is selected from the rising phase of the Solar Cycle, so as to include magnetic field structures 
representative of both Solar Maximum (in the newly-emerged active regions) and Solar Minimum 
(in the remnant regions and polar fields at higher latitudes). This period has the further advantage 
that it was used in our previous simulations to investigate the paramete r dependence of the mode l 
(most recently with regard to the formation and ejection of fiux ropes, lYeates &: Mackayll2009bl ). 
The start date of our comparison period is chosen to allow sufficient time for the sim ulation to 
evolve away from the initial condition on 1999 April 16 (see lYeates &: Mackavl l2009bl ) . while the 
end date just precedes the break in SOHO/EIT observations that occurred between 1999 September 
27 and 1999 October 5. 



^Available online from the CDAW data centre [http : // cdaw . gsf c . nasa . gov/ CME_list/ [ 
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In order to concentrate on well-observed CMEs, we ignored all events labelled "poor" by the 
CDAW observers. After studying the LASC0/C2 running difference movies of each remaining 
event, a further 12 very weak events were removed, as were 19 events which we could not clearly 
identify to be independent from other events. In addition, our study of corresponding EIT obser- 
vations (to be described below) led us to split one event in the CDAW catalog — on 1999 August 
18 at 05:54 — into two separate eruptions. A list of 330 CMEs remained, including both halo and 
limb events. 

There are some well-known uncertainties in the CDAW observations. Firstly, the y are sensitive 



to pr ojection effects, with CMEs in the plane of the sky being better observed (jHudson et al 
20061 ). However, since our data cover several solar rotations, there should be no systematic bias 



in the overall longitude distribution. Secondly, the selection of events is subjective. Recent ly, an 



automated CME catalog, CACTus (jRobbrecht &: Berghmans 



2004 



Robbrecht et al.l l2009al ) , has 



been developed to reduce subjectivity in the detection of CMEs by detecting radial motion in 
height-time maps of LASCO data. Although many additional events are detected in CACTus, 
they tend to be na rrow events, background outflows, or multiple detections of the same CME 
(jYashiro et al.l 120081 ). Thus the manual CDAW catalog is sufficient for comparison with simulated 
flux rope ejections. The third main limitation of LASCO observations is that they do not show the 
initiation locations of eruptions in the low corona, beca.use CMEs often rn ove non-radially before 



they reach the C2 field of view (e.q. JPlunkett et al. 



2001 



Attrill et al.ll2009l ). In this study we have 



used additional observations in extreme ultraviolet (EUV) to determine, where possible, the CME 
source locations in the low corona. Our unsuccessful attempt to use an alternative Ha data set is 
described in Section 12.21 



2.1. EUV Images 



A number of low-coronal features associated wi th CMEs appear in the E UV images from 
SOHO/EIT (Extreme ultraviolet Imaging Telescope; iDelaboudiniere et al.lll995l ). These include 



Table 1. Numbers of CMEs with associated EUV source locations. 



EIT source category 



Number of CMEs 



1 (clear) 

2 (plausible) 

3 (far-side) 

4 (no associated EIT source) 
No EIT data 

Total 



98 
44 
55 

124 
9_ 

330 
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Fig. 1. — Source locations associated with observed CMEs, in (a) latitude- longitude and (b) 
latitude-time. Asterisks show locations determined from (category 1) EUV observations, and dia- 
monds show locations determined from Ha data (Section 12. 2p . In (b) light-gray shading denotes 
LASCO data gaps recorded on the CDAW website. 



erupting filaments, f lares, post-eruption a rcades, and transient coronal dimmings, in addition to 



EIT coronal waves (jPlunkett et al.ll200ll ). BIT observations at up to 12-minute cadence were 
available in the 195A filter for the hours leading up to and following 321 of the 330 CMEs in our 
study period. Each of these were studied manually, in conjunction with the LASCO movie showing 
the CME. Possible CME signatures were recorded. We were unable to use EIT images in the 17lA, 
284A, or 304A filters as images were available only on one or two days during each month. To 
minimize subjective bias, two of us (ARY and GDRA) carried out this analysis independently, 
before comparing results and compiling a final event list (given in Appendix [A]). Further, each 
CME was assigned to one of the following categories: 



1. Clearly associated front-side source visible in EIT. 

2. Plausibly associated front-side source visible in EIT. 

3. Source becomes visible in EIT above the limb, but originated behind limb. 

4. No plausibly-associated EIT source evident. 



The first two categories give two degrees of certainty to our front-side source identifications. The 
third category describes CMEs where the low coronal source of the CME becomes visible in EIT 
above the limb, but its source on the solar disk lies behind the limb, so that a longitude position 
can not be determined. The locations of the possible sources in categories 1 and 2 were recorded 
manually by overlaying a latitude-longitude grid on the EIT movies. Identified source locations in 
category 1 are shown by asterisks in Figure [TJ Table [1] lists the number of CMEs in each category. 
For the comparison with simulations in this paper, we use only the 98 source locations in category 
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1, i.e., those clearly identified with CMEs. This represents only 30% of the observed CMEs, 
highlighting the difficulties associated with identifying CME source regions in the low corona. 
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Fig. 2. — Histograms showing (a) latitude, (b) longitude, and (c) time distributions of observed 
CME sources (gray shaded) and simulated ejections in run A4 (solid lines). Dashed lines show 
distributions of apparent latitudes and times in the original LASCO/CDAW data. Bin sizes are 5° 
in (a), 15° in (b), and 5 days in (c). 

Latitude, longitude and time distributions of CME source regions identified in EUV are shown 
by the filled gray histograms in Figure [2l The bir nodal latitude distribut i on of EUV sources is 
consistent with those found over the same period bv lCremades fc Bothmeii (|2004l . see their Figure 
10b), who carried out a similar SOHO/EIT analysis for a subset of "structured" CMEs between 
1996 and 2002. Although ou r analysis in this paper is restricted to the year 1999, the results of 
Cremades &: Bothmerl (j2004l ) indicate that the latitude distribution changes over the solar cycle, 
with an extension of the CME source distribution right down to the equator at solar maximum, 
in addition to a high latitude branch corresponding to polar crown filament eru ptions. This cycle 
variation is supported by observations of the latitudes of disappearing filaments (jPoioga &: Huang 
2OO3I ). 



The longitude distribution of CME sources is le ss structured, although from Figure Wl^a) we 
see a strong clustering into several activity complexes (jOaizauskas et al.lll983l ) , such as at longitude 
250° in the southern hemisphere, or 330° in the northern hemisphere. These active longitudes are 
also responsible for major fluctuations in the time distribution (Figure [2}:) . 

The dashed histograms in Figures [2)||a) and (c) show the apparent latitude and time distribu- 
tions of all 330 CMEs fror n the CDAW catalogue. A s expected, the distribution differs from those 
of the low coronal sources ([Gopalswamv et al.ll2003l ). The two main differences are (1) more events 
at t he equator, consistent with the deflection of some events towards the equator as they propagate 



out (Plunkett et al. 



2001 



Attrill et al.ll2009l ). and (2) more events at high latitudes in the Northern 



hemisphere. This north-south asymmet ry in the latitude dis tribution of LASCO CMEs is a feature 
of the year 1999, evident in Figure 3 of lYashiro et al.l ([200J). 
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The solid black histograms in Figure [2] show the equivalent distributions of flux rope ejections 
in simulation run A4. The simulations will be described in Section [3l 



2.2. Ha Event Database 



As an alternative to using EUV images, we als o tried identify i ng sou rce locations of observed 
CMEs using a database of Ha events assembled by Howard et al. (j2008l ). This database lists the 



locations of flares and disappearing fllaments in Solar Geophysical Data that may be associated 
with CMEs in the CDAW catalogue. A surface event is associated with a particular CME if it 
occurred within ±1 hour of the CME onset time (estimated assuming a constant outward speed), 
and if its latitude and longitude are in the same quadrant as any part of the angular span of the 
CME. It should be noted that, in this database, some CMEs are associated with multiple possible 
surface events, and many with none. Also, some surface events are listed with multiple CMEs. In 
an attempt to select only those CMEs with reasonably consistent locations given by the Ha events, 
we filter out CMEs whose list of associated Ha events has standard deviations greater than 10° in 
latitude or 20° in longitude. This leaves 137 CMEs with consistent locations in the Ha database, 
in the sense that the Ha events associated with each CME are reasonably close together spatially. 
These source locations are shown by diamonds in Figure [H It is clear from Figure [U that, like the 
locations identified in EUV, those identified in Ha data are not uniformly distributed, but rather 
are clustered into several major activity complexes. 

Unfortunately, we find that the Ha source locations are not consistent with the EUV source 
locations. Firstly, it is evident from Figure [T] that the Ha sources are even more strongly clustered 
into active regions, with no sources whatsoever observed at latitudes above about 40°. The corre- 
lation coefficient between the EUV and Ha spatial distributions is only 0.68 (when binned in 20° 
latitude bins and 30° longitude bins, following the method described in Section H]). The discrepancy 
becomes even more serious when considering the source locations identified with particular CMEs. 
There are 42 CMEs with both a category 1 source location in EUV and a location identified in 
Ha. However, the Ha location overlaps the EUV location (to within ±10° in latitude and ±20° in 
longitude) in only 22 of these cases. We are inclined to favor the EUV locations rather than the Ha 
locations (at least for category 1 EUV events) because the EUV data were examined in detail on a 
case-by-case basis, which differs from the approach used to compile the Ha listing. The less careful 
association of Ha events with individual CMEs is exemplified by 23 events which were clearly seen 
to originate from behind the limb in EUV, yet in the Ha database were associated with flares that 
happened to occur co-temporally on the solar disk. The Ha observations do have the potential 
advantage of a higher cadence than the EUV images (minutes rather than over ~ 10 minutes), so 
could in principle identify additional CME sources that were missed in EUV. However, in this par- 
ticular database, we have no clear way of selecting which of many possible Ha events are actually 
associated with each (or any) CME. Therefore, we base the comparison with our simulations on 
the EUV data only. 
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3. Coronal Magnetic Field Simulations 



Ou r numerical simula t ions of the large-scale coronal magnetic field evol ution were described in 



detail in lYeates &: Mackavl (|2009bl ) . and are based on the mean- field model of lvan Ballegooijen et al 
([20001) • Briefly, the large-scale mean magnetic field Bq = V x Aq evolves via the non-ideal induction 
equation 

voxBo-r/jo, (1) 



dAo 



dt 



where the turbulent diffusi yitv H is given by a background va lue t/q and an enhancement in regions 
of strong current density (IMackav &: van BallegooiienI l2006l ) . The velocity is determined by the 
magneto- frictional technique (jYang et al.lll986l ) as 

1 jo X Bo 



vo = -- 



+ ^out(^)r, 



(2) 



where the first term approximates the evolution as a sequence of quasi-static force-free equilibria, in 
response to flux emergence and shearing by large-scale surface motions on the lower, phot ospheric, 
boun dary. The surface motions are given by the standard surface flux transport model (jSheelev 
20051 ). The second term in equation ([2]) represents a radial outflow imposed only near the upper 
boundary (r = 2. Si?©), to represent the effect of the solar wind in radially opening magnetic field 
lines. 



In this study we use the simulation runs described in lYeates &: Mackavl (|2009bl ). All cover 
the same time period, starting on 1999 April 16. The initial condition is a potential-field source- 
surface extrapolation taken from a synoptic normal-component magnetograin froin the National 
Solar Observatory, Kitt Peak, corrected for differential rotation (jYeates et al.l 120071 ). The surface 
and coronal magnetic fields are then evolved continuously for 177-days. During this evolution, 
119 new bipolar magnetic regions are inserted into the simulation, with location, size, tilt angle 
and magnetic flux determined from Kitt Peak synoptic maen etograms ([Yeates et al.l 120071 ) . The 
bipolar regions take the idealized mathematical form given in lYeates et al.l ( 2008a ). Because the 
simulation is non-potential, non-zero currents and magnetic helicity are generated in the corona 
during the evolution. This arises not only due to shearing by the large-scale surface motions, but 
also because the newly-emerging bipolar regio ns may be given a, non-zer o helicity, controlled by 
a "twist" parameter /?, described in detail by lYeates k. Mackavl (j2009bl ). As discussed in that 
pap er, although techni ques to measure the twist in observed active regions have been developed 
(see iNandv et al.l |2008| ) , available observations are limited such that we cannot reliably determine 
the amount of helicity and thus optimum value of f3 to model each individual active region. For 
simplicity, we have assumed that all bipoles in each hemisphere have the same value of /3, but with 
different val ues of in each hemi sphere (taking into account the hemispheric dependence of helicity 
observed by iPevtsov et al.l Il995l ) . We have run a series of simulations with different values of (5 
to study the effect of this parameter on the formation and evolution o f magnetic flux ropes. Th e 
different runs are summarized in Table [21 repeated for convenience from [Yeates &: Mackavl (|2009bl ) . 



A natural consequence of the quasi-static evolution simulated in this model is the accumulation 
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above polarity inversion lines of a xial magnetic field, in the form of twisted magnetic flux ropes 



van Ballegooiien &: Martenslll989l ). The formation o f flux ropes in the mean-field mod el has been 



studied in detail for a simple 2-bipole configuration (IMackav &: van BallegooiienI 120061 ) . and more 
recently in the global simulations used in this paper (jYeates &: Mackavll2009bl ). In these simulations 
it is found that, once the axial magnetic field in a flux rope grows too strong relative to the 
overlying field, the flux rope loses equilibrium, rises, and is ejected through the top boundary of the 
computational domain. It is these flux rope "ejections" that we compare in this paper to observed 
low-coronal CME source regions. 



Yeates &: Mackavl (j2009bl ) developed automated methods to identify flux rope structures and 
their ejection in the global simulations, allowing objective comparison between different simulation 
runs. In Figure [3] we show the results of this procedure for one particular day (1999 August 28, day 
of year 240), in simulation run A4. Figure [3l[|b) shows the simulated magnetic field viewed from 
Carrington longitude 130°, with the radial magnetic field on the solar surface shown in grayscale and 
coronal field lines in blue (if closed) or orange (if open). For comparison, the Kitt Peak magnetogram 
for CR1953 is shown in Figure [3]^a). Most field lines in Figure [3lfb) have been tr aced from the fiux 
rope points selected by the automated technique (see lYeates Sz Mackavl l2009bl ) . These flux rope 
points are shown projected on the solar surface in Figure [3l^c) , where each flux rope structure is 
numbered. The background shading again shows the radial surface magnetic fleld on the same day 
(here white is positive and gray negative). Figure [3)||d) superimposes all flux rope points selected 
between day 235 and day 245, with the points colored red if they are involved in an ejection during 
this period, blue otherwise. 

The number of ejections in the simulations is found to depend (to some extent) on the twist 
of emerging bipolar regions, with more ejections for /? of greater magnitude, or with the observed 
majority sign in each hemisphere. The right-most column of Table [2] shows the rate of flux rope 
ejections in each simulation run, over the period selected for the observational comparison in this 
paper (1999 May 13 to 1999 September 26). In the next section, we compare the distribution of 
these flux rope ejections with observed CME source locations. 



4. Comparison between Simulations and Observations 

Having identifled the locations both of a subset of observed CME sources (Section [2]) and of 
simulated flux rope ejections (Section [3]), we now compare the distributions of each in latitude, 
longitude and time. Of course, given the limited number of observed CMEs where sources could 
be reliably determined, this comparison is necessarily approximate. 

Figure [4] shows the locations of observed CMEs and simulated ejections (run A4), for each 
Carrington rotation during the simulation period. Here the crosses show locations and extents of 
EIT category 1 CME source locations, while clusters of different colored points refer to different flux 
rope ejections in the simulation. Each color is equivalent to a separate cluster of red points in Figure 
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Fig. 3. — Overview of the simulated magnetic field and flux rope detection: (a) normal-component 
magnetogram from NSO/Kitt Peak for CR1953; (b) snapshot of the simulated magnetic field on 
day 249; (c) locations of magnetic flux ropes detected on day 240; (d) locations of flux rope points 
(blue) and ejections (red) detected between days 235 and 245. 



[3l^d). Note that we only compare simulated ejections falling within the comparison period, which 
begins midway through CR1949 and ends midway through CR1954, hence the smaller number of 
observed CMEs shown in Figures Hl^a) and (f). Figure H] shows that the simulated ejections clearly 
do not match the observed locations on a one-to-one basis. This is not surprising, given both the 
limited observations and poorly constrained simulation parameters such as the twist of emerging 
bipolar regions, their date of emergence, or the turbulent diffusivity in the corona. However, 
the locations of many simulated and observed ejections do overlap, and there is a region in the 
Southern hemisphere around Carrington longitudes 150° with few ejections in either observations 
or simulations. Furthermore, the overall latitude distribution for run A4 is bimodal and broadly 
consistent to that of the observed source locations (Figure [2^) . 

To quantify the association between observed and simulated distributions, we carry out a 
straightforward correlation analysis. Each list of ejection locations (simulated or observed) is binned 
in latitude, longitude and time. The bin sizes in latitude and longitude are chosen so as to take 
into account the spatially-extended nature of the source regions, which are not single points. From 
cumulative distributions of the latitudinal and longitudinal extents of simulated ejections and ob- 
served source regions, we select bin sizes of 20° in latitude and 30° in longitude, so as to be larger 
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Fig. 4. — Carrington maps showing locations of simulated ejections in run A4 (clusters of colored 
points), and category 1 CME source locations observed in EUV (with crosses showing the approx- 
imate extent of the source region). The polarity of the photospheric radial magnetic field in the 
simulation is shown in the background (white for positive, gray for negative), for the middle day 
of each Carrington rotation. 



than 80% of events. The variation of bin size in time will be considered below. 

To illustrate the technique of estimating a quantitative correlation, the binned distribution of 
source locations for simulation run A4, in latitude and longitude, is shown in Figure [5^a). Here 
black indicates no events in that bin (over the whole 136 days), with white indicating the most 
events. Figure [5{b) shows the equivalent distribution for the observed CME sources. To assign a 
quantitative correlation, we compare the number of simulated ejections with the number of EUV 
sources in each bin. This is shown by the scatterplot in Figure [5|^c) , where the size of each circle 
indicates the number of bins with those numbers of simulated and observed sources. The Pearson 
linear correlation coefficient is then computed. In this case it is 0.49, indicating a significant, but 
not particularly strong, positive correlation. 
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Fig. 5. — Example analysis of correlation in spatial distributions. Histogram in (a) shows the 
distribution of simulated ejections (run A4), and (b) shows the distribution of observed EUV 
sources. The correlation between the number of events in each bin is shown in (c), where the 
sizes of the circles are proportional to the number of bins containing those numbers. The linear 
correlation coefficient is 0.49. 



4.1. Spatial distributions 

Table [3] shows the linear correlation coefficients between the binned latitude-longitude distri- 
butions of flux rope ejections in the various simulation runs and that of the observed CME sources. 
All simulation runs show a positive correlation except for run AN (with no emerging bipoles). The 
correlations are stronger for the simulation runs with the observed majority sign of emerging bipole 
twist in each hemisphere (negative in the Northern hemisphere and positive in the South, i.e., runs 
A2, A4, A6 and V4) than for the runs with either untwisted bipoles (AO) or bipoles with the oppo- 
site sign of twist (runs Am2, Am4 and Am6). However, this is a secondary effect and the positive 
correlation is present in all runs except AN (where there are no emerging bipoles). The highest 
correlation coefficient attained is 0.49 for run A4. Interestingly, this is the run that best agreed 
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Table 2. Summary of different simulation runs. 



Run 


/3 in N. hemisphere 13 in 


S. hemisphere 


770 (km^ s ^) 


Ho (kms^^)'' 


Flux rope ejections per day'' 


AN 


No emerging reg 


ions 


45 


100 


0.17 ±0.03 


Am6 


0.6 


-0.6 


45 


100 


1.09 ±0.16 


Am4 


0.4 


-0.4 


45 


100 


1.02 ±0.15 


Am2 


0.2 


-0.2 


45 


100 


0.64 ±0.10 


AO 








45 


100 


0.72 ± 0.11 


A2 


-0.2 


0.2 


45 


100 


0.99 ±0.15 


A4 


-0.4 


0.4 


45 


100 


1.15 ±0.17 


A6 


-0.6 


0.6 


45 


100 


1.27 ±0.19 


D4 


-0.4 


0.4 


22.5 


100 


1.46 ±0.22 


V4 


-0.4 


0.4 


45 


50 


1.12 ± 0.17 



^Peak value of radial outflow velocity, at r = 2.5-Rq. 

''Number of ejections per day between 1999 May 13 (da y of year 133) and 1999 S eptember 26 (day of year 269). 
Errors are those for the automated flux rope detection (see Yeates fc Mackav 2009bl ) . 



Table 3. Correlation between spatial distributions of observed CME sources and simulated 

ejections. 



Run Correlation coefficient 



AN 


0.08 


Am6 


0.37 


Am4 


0.33 


Am2 


0.33 


AO 


0.33 


A2 


0.42 


A4 


0.49 


A6 


0.43 


D4 


0.31 


V4 


0.40 
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with observations of filament chirality in our earlier comparison (jYeates et al.ll2008al ). Although a 
correlation of 0.49 may seem low, it is still significant, given the large-scale simplified nature of the 
simulations and the uncertainties associated with the observations. 



4.2. Correlation in space and time 




Time Bin Size (days) 



Fig. 6. — Linear correlation coefficients between simulated ejection locations and observed CME 
source locations as a function of bin size in time. Asterisks/solid line show simulation run A4, 
diamonds/dashed line show run AN, and dotted lines show other simulation runs. 

Given that the overall latitude-longitude distributions of simulated flux rope ejections and 
observed CME sources are positively correlated (though not very strongly), are these correlations 
maintained if the data are also binned in time? Figure [6] shows how the correlation coefficients 
decrease as the bin size in time is reduced from the whole comparison period (136 days) — which 
corresponds to the spatial correlations in Table [3] — to shorter and shorter bins. We see that signifi- 
cant positive correlations between the simulated ejections and the observed sources are found only 
with bin sizes of 34 days or longer. This is readily understandable because the observational in- 
put driving the simulations — synoptic magnetogram data — is only updated every 27 days. Clearly, 
if the model is to reproduce observations on a shorter timescale then more frequent updating of 
emerging flux would be required. 

Note that the timescale of flux rope formation and loss of equilibrium in our model also depends 
on the turbulent diffusivity in the corona, which is not directly constrained by observations. In run 
D4 we halve the value of the coronal diffusivity, which results in larger flux ropes and 25% more 



ejections. However, as pointed out in lYeates &: Mackavl (|2009bl ). many of the flux ropes in run D4 



are rather highly twisted to be realistic. Furthermore, Table [3] shows that the spatial distribution 
of simulated ejections in run D4 is less well correlated with observed CME source locations (with 
a correlation coefficient of 0.31, compared with 0.49 for run A4). 
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5. Discussion 

Having identified some (limited) correlation between the simulated flux rope ejections and 
observed CME sources on timescales of a month or longer, we now consider the origin of this 
correlation. 

The locations where flux ropes form in the simulation are determined by the structure of the 
magnetic field. F lux ropes form above pol arity inversion lines where axial fields build up after fiux 



cancellation (see lYeates Sz Mackavi l2009bl ) . The axial (sheared) components originate both from 
the emergence of twisted bipoles and from shearing by large-scale surface motions, so that the 
distribution of flux ropes in the simulations depends to some extent on the balance of these two 
contributions. 

From the dashed line in Figure El it is clear that the correlation disappears in simulation 
run AN, where no new bipolar regions were inserted after the initial condition. This suggests 
that the correlation observed in the other simulation runs arises primarily from the magnetic 
field distribution on the solar surface. This is essentially captured by the surface flux transport 
component of the simulations. To investigate this idea further, we consider purely the distribution 
of Br on the solar surface r = Rq, which is the same in all simulation runs except AN. This 
distribution represents the observational input driving the coronal magnetic fleld evolution in our 
model. It is summarized in Figure [T^a), which shows the latitude- longitude distribution of {\Br\), 
where the average is taken over all 136 days in the comparison period. The magnetic field is 
clearly non-uniform over the solar surface, and is concentrated in several major activity complexes, 
corresponding well to the clusters of observed CME sources identified in EUV (shown by asterisks) . 
For comparison, the triangles show the locations of flux rope ejections in run A4, which are more 
evenly spread over the solar surface than the observed source regions. In a similar vein. Figure 
[Tljb) shows a map of {\'VhBr\), the time-averaged horizontal gradient of Br on the solar surface. 
The distribution resembles that of {\Br\), except that there is a greater concentration in the active 
region belts relative to higher latitudes. 

Table 4. Correlation between observed CME sources and simulated surface magnetic field. 



Time bin size (days) 


Correlation with ( -Br ) 


Correlation with ( V^-Br ) 


136 (full period) 


0.62 


0.67 


68 


0.50 


0.53 


34 


0.41 


0.43 


17 


0.30 


0.33 


8 


0.23 


0.24 


4 


0.18 


0.19 
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Fig. 7. — Locations of simulated ejections and observed CME sources, overlaid on plots of magnetic 
field properties averaged over the comparison period: (a) radial surface field strength {\Br\); (b) 
horizontal gradient in radial surface field {{S/fiBrl); (c) current helicity (|a|) at height r = I.OSRq 
in run A4. In each case asterisks show observed CME sources and triangles show locations of 
simulated ejections in run A4. 

We may consider using these quantities {\Br\) and {\VhBr\) themselves as predictors of CME 
source locations. For this purpose we use the same binning as before, calculating the mean values 
of these new quantities in each bin. Table [4] shows the resulting correlations with observed CME 
source locations, for various choices of time bin size. Notably, both {\Br\) and {\VhBr\) are better 
correlated with observed CME sources than are the flux rope ejections in any of the simulation 
runs. For example, the spatial distribution (136-day bins) of EUV sources has a correlation coeffi- 
cient of 0.62 with (I-Bri) and 0.67 with (|V/,i3r|), but only 0.49 with flux rope ejections in run A4. 
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From Table IH we see that the correlations are maintained for time bins as short as 17 days, and 
are consistently better than the simulated ejections (Table [3]). Table H] also shows that {\VhBr\) 
is better correlated with observed CME locations than is (|-Br|)- This is presumably because the 
former is more concentrated in active regions, where the majority of observed CME sources are 
located. Indeed, measures of active region complexity based on strong gradients of radial pho- 
tospheric m agnetic flux ha ve been found to correlate with pro duction from those regions both of 
large flares (jSchrijveiil2009l ) and of CMEs ([Falconer et alJl2008l ). If vector magnetograi ns are avail- 



able, a lternative characterisations of the magnetic complexity are possible; for example, iHahn et al 



(j2005l ) find a correlation between magnetic twist and the locations of Ha flare signatures within 
active regions. Of course, although {\VhBr\) is better correlated with the distribution of observed 
CMEs, it could not be used to predict their time of occurrence, because the photospheric magnetic 
field does not change significantly between the pre-eruption and post-eruption state. Therefore, 
time-dependent models of the 3D coronal magnetic field evolution must still be developed for this 
purpose. 



Yeates &: Mackay (j2009bl ) showed that the amount of helicity in emerging bipolar regions 



does alter the simulated flux rope ejections, but this would appear to have only a small effect on 
the correlation with observed CME sources. Figure d^c) shows the time-averaged cu rrent helicity 
a = j • B/i?^ in simulation run A4, at a height of r = l.OSi?© in the low corona (see lYeates et al 



2008bl . for a discussion of the distribution of a). From the triangles in Figure [Tl^c), we see that 
the ejections in run A4 always occur at locations of strong in the simulation. However, if we 
compute the correlation between the binned distribution of (|a|) and the observed CME sources, 
we find a correlation coefficient of only 0.38, much lower than for (|-Br|) or {\VhBr\), and also 
lower than for the simulated ejections. This is because the observed CME sources are much more 
strongly concentrated in active regions than is the distribution of current helicity in the simulation. 
However, there are some EUV sources observed at higher latitudes, mainly quiescent filament 
eruptions. So the clustering of observed CME sources at active latitudes does not preclude the 
presence of coronal currents at higher latitudes. It is possible that there are two essentially different 
subsets of CMEs: those produced by quasi-static shearing and flux cancellation, not dissimilar to 
the distribution of flux rope ejections in our model, and a second set of events concentrated in active 
regions and resulting from more dynamic and/or smaller-scale processes not included in our model. 
Indeed, observations over several decades suggest that CMEs originating from active regions and 
associated with flares have different kinematic properties from those associated with prominence 
eruptions outside active regions. Th e former are faster, while the latter start more slowly and 



show gradual acceleration with height (jCosling et al 



1976 



MacQueen Sz Fisheiill983l : ISheeley et al 



1993). However, it is still not c ertain whether CMEs with different speeds are qualitatively different 
(jCopalswamv et~aL 2006 ). and Low Sz Zhang ( 2002 ) show how both could arise from flux ropes in 
different configurations. The EUV observations in this paper do not generally allow us to directly 
establish the presence or absence of a flux rope morphology, except for well-observed filament 
eruptions. 
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6. Conclusion 

In this paper, magnetic flux rope ejections in a global model of the coronal magnetic field 
evolution have been compared directly with low-coronal CME source regions identified in EUV, over 
a 136-day period in 1999. Despite limitations in identifying the source regions of observed CMEs, 
which was unambiguously possible for only a third of events, we find some definite correlation 
between the distributions of simulated ejections and observed CME sources. However, there are 
also clear differences between the observed and simulated distributions; the model produces fewer 
ejections overall, and much less clustering of ejections in active regions, even given the incomplete 
observations. 

The correlation that is present between simulated and observed ejections originates primarily 
in the large-scale distribution of magnetic field on the solar surface, the basic observational input 
driving the coronal simulations. The correlation is only significant when the ejections are binned 
in time bins of at least 34 days, which is approximately the frequency (27 days) at which new 
observations of the surface magnetic field are are available as input to coronal simulations. There 
is a weak dependence of the correlation on the helicity in emerging active regions, but this is a 
secondary effect. Interestingly, the choice of emerging helicity that performs best is simulation run 
A4, which is also the run that most accurately reproduced the observed pattern of filament chirality 



(jYeates et al.ll2008al ). However, a better predictor of CMEs than the locations of fiux rope ejections 
is found to be simply the radial magnetic field on the solar surface, and in particular its horizontal 
gradient. This quantity is readily computed from surface flux transport simulations alone (or, 
periodically, from observed line-of-sight magnetic maps). Its good correlation with the observed 
CME source locations arises because the observed sources are clustered strongly into active regions, 
where the magnetic field is strong. By contrast, the simulated flux rope ejections are more evenly 
distributed over the solar surface, following the distribution of current helicity in the model. They 
do show some concentration towards active latitudes, reflecting the smaller scales of current helicity 
found in the more complex magnetic fields there, but the distribution is clearly different. 

Our results allow us to place constraints on the possible initiation mechanisms of CMEs, 
at least in the rising phase of the 11-year Solar Cycle. Since flux ropes in the simulations are 
formed by gradual, quasi-static shearing of the magnetic fleld by large-scale motions, along with 
flux cancellation, we conclude that this mechanism cann ot be responsible for all CMEs o n the Sun, 



although it is sufficient to produce a signiflcant fraction (jLowll200ll : iGibson &: Fanll2006l ). Whether 
this fraction varies over the Solar Cycle is a question for future research. Based on our results 
here, we propose that there are two subsets of CMEs produced by essentially different mechanisms: 
those produced by large-scale, gradual transport of helicity (which may originate either from active 
regions or shearing by surface motions), and a second population concentrated in active regions. 
The former occur on timescales of weeks or even months, and are simulated in our model. The 
latter must occur on much shorter timescales, connected with the dynamic emergence of magnetic 
flelds in active regions and their energetic restructuring. Since our model was originally developed 
to study the large-scale transport of magnetic helicity in the corona, it cannot produce the second 
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type of CME using the present form of input data (synoptic magnetograms) . To reproduce the 
frequent ejections observed from active regions wih require, as a minimum, magnetogram data at 
much higher spatial and time resolution for input to the quasi-static simulations, at least within 
active regions themselves. However, given the dynamic nature of many events, it is likely that 
detailed, time-dependent MHD simulations of complex struct ures will ne ed to be performed, based 
on observations of flux emergence such as those by Hinode (ILited l2009l ) . At present, developing 
a detailed understanding of the initiation of such CMEs that could possibly lead to predictive 
capability is far out of reach, although there have been recent advances in modelling flux emergence 
in simplified configurations (| ArchontisI l2008l ) . 

Perhaps a more achievable goal in the short term is to understand the evolution leading to 
the initiation of the first type of CME, those produced by quasi-static shearing. Above the active 
latitudes, remnant magnetic fields from active regions do not retain most of the complexity of 
their originating regions, and shearing motions along with flux cancellation over longer time-scales 
become more important. Our existing model provides a starting point to understand the net 
effect of these motions on flux rope development. However, the transport of helicity over the solar 
surface means that, at least at mid- latitudes, the formation of flux ropes is sensitive to the amount 
of magnetic helicity in individual emergin g active regions, which is not yet routinely observed, 
despite recent improvements in techniques (jPemoulin &: PariatI l2009l ) . The forthcoming launch of 
NASA's Solar Dynamics Observatory (SDO) should improve the situation, with full-disk vector 
magnetograms at high cadence and resolution. This should allow us to put better constraints on 
the amount of helicity in the solar corona, and to further validate our quasi-static model for flux 
rope formation. There is also an open question as to the origin of sheared fields at high latitudes. 
In our present simulations, where helicity emerges from the solar interior only in active regions, 
differential rotation is the dominant so urce of helicity at h igh latitudes, but produces the opposite 
sign of helicity there to that observed (jYeates et al.ll2008bl ). 

Finally, we note that the observational identification of CME source regions is not straight- 
forward. Our attempt to use an alternative database of CME source regions in Ha highlights the 
need for careful case-by-case investigation, preferably in multiple wavelengths. With a maximum 
cadence of 12 minutes in EIT 195A data, we were able to confidently determine the CME source 
regions of only 98 events on-disk, with a further 55 originating from behind the limb, out of a 
total of 330 CMEs. This is certainly incomplete. The situation should be improved in future, 
using higher-cadence observations at multiple EUV wavelengths, for example from SDO. However, 
if future aims are to identify the relative importance of different CME initiation mechanisms in the 
global context, then a comprehensive catalog of CME sources must be built up over a long period. 
Only then can more definitive comparisons with theoretical models be made. Again, SDO offers 
some promise in this respect: for example, automated feature- finding algorithms fo r detecting the 



low-coronal signatures of CMEs in the anticipated SDO data are in preparation (jMartens et al 



2009; 



Attrih fc Wills-Davevll2009l ;i. 
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A. Observational Event List 

Table [5] gives the list of 330 LASCO CMEs, with associated EUV source information where 
a source was found. This includes source quality, approximate location, and notes on the type of 
signature(s) in each case. 
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Table 5. Low coronal source regions for LASCO CMEs. 



LASCO CME Data'' Low Coronal EUV Source 



Date 




Time 


CPA'' 


Width'^ 


Quality'* 




itude Latitude 


Notcs° 
















Min 


Max Min 


Max 










(UT) 


fdeff) 


V &/ 












1999 


May- 


13 


23:26:07 


311 


105 


2 


50 


70 15 


25 


PEA near NW limb 21:35 


1999 


May 


14 


03:06:05 


55 


63 


2 


-90 


-70 40 


50 


Activity & dimming over NE limb 02:58 


1999 


May 


14 


06:50:05 


234 


46 


1 


40 


90 -50 


-30 


Large QF eruption 03:22-f 


1999 


May 


16 


17:51:35 


33 


45 


2 


60 


90 -20 


10 


Activity in large loop structures on limb, 
maybe partial F eruption 


1999 


May 


17 


00:50:07 


293 


113 


3 




far-side NW 




Eruption-opened loop structures above limb 


1999 


May 


19 


23:02:31 


274 


17 


4 










1999 


May 


19 


23:02:31 


57 


94 


4 










1999 


May 


20 


08:26:05 


47 


69 


4 










1999 


May 


20 


13:28:21 


313 


119 


1 


-20 


25 


50 


Clear eruption, dimming, & PEA loops from 
centre NH, 08:48-1- 


1999 


May 


20 


16:26:05 


21 


111 


1 


-45 


30 


65 


Restructuring in NH AR after previous euption, 
PEA from 18:35 


1999 


May 


20 


21:26:08 


115 


91 


4 










1999 


May 


21 


10:50:05 


341 


130 


4 










1999 


May 


21 


17:50:06 


20 


218 


3 




far-side NE 




PEA above NE limb 20:11 (also clear in He304) 


1999 


May 


23 


07:40:05 


288 


67 


? 










1999 


May 


23 


19:06:01 


50 


21 


7 










1999 


May 


24 


10:30:05 


101 


28 


7 










1999 


May 


24 


17:07:31 


39 


44 


4 










1999 


May 


25 


05:06:05 


81 


12 


4 










1999 


May 


25 


07:26:51 


103 


35 


3 




far-side SE 




PEA above limb 07:35-f 


1999 


May 


25 


10:50:05 


268 


178 


1 


70 


90 -50 


-10 


Clear eruption on SW limb 10:23, dimming, 
coronal wave, flare 


1999 


May 


25 


16:26:05 


26 


133 


4 










1999 


May 


25 


23:26:05 


5 


138 


4 










1999 


May 


26 


04:26:05 


44 


71 


2 


-45 


-35 20 


30 


Surge-like small eruption -|- flare from NE sector 02:35 


1999 


May 


26 


05:26:05 


106 


51 


3 




far-side SE 




Dimming & PEA over SE limb 05:11+ 


1999 


May 


26 


08:06:05 


321 


101 


1 


35 


90 40 


70 


Slow QF eruption near NW limb starts 23:47 on 
May 25, dynamic phase starts 07:24, dimming, PEA 


1999 


May 


26 


20:26:05 


38 


17 


2 


-60 


-40 15 


30 


Strong brightening in AR near NE limb 19:43-1- 


1999 


May 


26 


22:26:05 


110 


70 


1 


-90 


-35 -60 


-20 


QF on SE limb starts to erupt 19:14, 



slow eruption until disappearance at 22:35 



''As obtained from the CDAW catalogue at |http : //cdaw.gsf c .nasa. gov/CHE_list/[ 

^'Central position angle. 
'^Angular width in plane of sky. 

''Quality 1-4 of the EUV source identification, as defined in Section[2] The symbol ? indicates that FIT data were unavailable 
at that time. 

''Abbreviations used include AR (active region), PEA (post-eruption arcade), and QF (quiescent filament). 

Note. — Table |5] is published in its entirety in the electronic edition of the Astrophysical Journal (or may obtained by email 
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request to the first author). A portion is shown here for guidance regarding its form and content. 
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